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- TITLE minesesent - 
IDENT /1-0 ; File: MTHCDSQRT.MAR Edit: SBL1003 


MASRBAAAAAASAASAASLALASALAAAAASALALALAALE SALES ASAA SAAS ARES ALES ALAR SESE E SERRA SES OS 


COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


o 
® 
& 
& 
® 
THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THiS SOFTWARE. OR OTHER * 
COPIES, THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY © 
OTHER PERSON. NO TITLE TO AND OUNERSHIP OF THE SOFTWARE 18 HEREBY 
TRANSFERRED. * 
® 
THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
CORPORATION. * 
® 
ee 
ty 
a 
® 
® 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


. . . . . . * o * . . * . 
* e@neeeneneneeneennnennennannne 


(AR RRASASAAASAASASAALALASASAREAALAALA SALAS ARRA ERR A RASA ESAS AAAS ASA AAAS ASAE SD OS 


“ABSTRACT: 
This module contains routine MTHS$CDSQRT - compute 
D COMPLEX*16 square root. 


VERSION: 1 
HISTORY: 


AUTHOR: 
Steven B. Lionel, 24-July-1979 


| 
(FACILITY: MATH LIBRARY | 
MODIFIED BY: 
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49 -SBTTL HISTORY ; Detailed Current Edit History 


; Edit History 


: 1 1-09) - Adapted from MTH$CSQRT version 1-003. SBL 24-July-1979 
; 1-002 = Change shared external reference to G* RNH go~sep-8t 
; 1-003 = Use general mode addressing. SBL 30-Nov-1981 


ee tt tN LL CC CC CC CC LL LLL LL 


o V04-00 


og | 
MTHSCDSQRT 16*SEP=1984 01:07:44 VAX/VMS Ma P | 
1-008 3Be 91:95:88 MTMRTL.SRCIMTHCDSQRT.MAR:1 9” “4 


DECLARATIONS 6-SEP=1984 11 


p 28 -SBTTL DECLARATIONS 
60 ; 
: 61 ; INCLUDE FILES: 
0 6¢ 3 
000 6 
000 64 ; 
000 65 ; EXTERNAL SYMBOLS: 
000 66 -DSABL GBL 
0000 67 -EXTRN MTHSCDABS 
0000 «= 68 .EXTRN MTHSDSQRT_RS 
$600 96 
0 0; 
0000 71 ; MACROS: 
0000 i 3 NONE 
0000 7 
0000 74; 
0000 75 ; PSECT DECLARATIONS: 
00000000 76 
0000 77 
0000 78 ; 
0000 79 ; EQUATED SYMBOLS: 
0000 80 
0000 81; 
0000 Hf 3: OWN STORAGE: 
0000 83 ; NONE 


| 
.PSECT _MTHSCODE PIC, SHR, LONG, EXE, NOWRT 


F 
MTHSCDSQRT 16-SEP-1984 01:07:44 VAX/VMS Macro v04-00 
1008 MTHSCDSQRT = compute D COMPLEX*16 square 6=-SEP-1984 91190588 MTHRTL.SRCJMTHCDSORT.MAR; 1 


5 -SBTTL MTHSCDSQRT = compute D COMPLEX*16 square root 


ae 
FUNCTIONAL DESCRIPTION: 


The square root of a complex number (r, i) is computed 
as follows: 


ROOT = SQRT((ABS(r) + CABS((r, 1))) / 2) 
Q= i / (2*ROOT) 


r j CSQRT((r, i)) 
>=0 an (ROOT, Q) 

<0 > (Q, ROOT) 

<0 <0 (-4, -ROOT) 


; CALLING SEQUENCE: 
CALL MTHSCDSQRT (result.wdc.r, arg.rdc.r) 
; INPUT PARAMETERS: 


arg = 8 ; The D COMPLEX*16 argument, passed 
; by reference. 
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00000008 

IMPLICIT INPUTS: 
NONE 

OUTPUT PARAMETERS: 


result = 4 3; The D COMPLEX*16 result, passed by 
3; reference. 
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00000004 0000 
; IMPLICIT OUTPUTS: 
NONE 


; COMPLETION CODES: 
NONE 


So 
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SIDE EFFECTS: 
SS$_ROPRAND If either part of argument is reserved operand. 


007C -ENTRY MTHSCDSQRT, “M<R2, R3, R4, RS, RO> 


MTHSFLAG_JACKET ; flag as math routine 


MOVAB G*MTHSSJACKET_HND, (FP) ¢ 
3; set handler address to jacket 
3; handler 
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52 08 BC 70 9 137 MOVD @arg(AP), R2 3 R2-R3B =r 
see Bee i A Cae 1A ees 
ar ; Put address of arg on stac 
ooovo000'Gr” OF FB 0015 140 CALLS. #12 C*MTHSCDABS i RO-RI = CABS((r, 9)9 
3 6 ic 141 ADDD R2, R RO s aed = ABS(r) + CABS(( i)) | 
0 QO 64 iF 196 MULD #0.5, RO 3 “R1 = (ABS(r) + CABS ((¢ 1} 7 2 
<r f- GF 1 0 § 14 JSB G*TASDSORT_RS ; ee = ROOT = SORT (aboves 
AC p 8p 144 MOVL arg(AP), R2 3; R2 => (r, i) | 
50 O2c 145 TSTD RO 3 je ROOT zero? 
06 12 00 146 BNEQ 18 : go ahead 
55 7C 00 147 CLRD R5 $ ate zero quotient 
98 11 88 138 BRB 2$ ; skip divide 
55 08 a 67 0034 149 1$: DIVD3 Re. 8¢Rg), R5 ; RS = i / ROOT 
5 00 $ 0039 150 MULD2 #0.5, R : RS =Q<2=i/ (2 * ROOT) 
8 3 003C 151 2$: STD (R2)+ :; if r positive 
15 18 00 138 BGEQ  RETRN : then return (ROOT, Q) 
3s ww 4 004 15 MOQ RO, R3 ; else switch oe and Q 
62 73 0043 154 TSTD = CR) : if i posit 
08 18 0045 #155 BGEQ RETRN1 3 then recure fa, ROOT) 
0 55 72 0047 156 MNEGD R5, R : else negate ROOT and Q 
5 53 72 QO4A 157 MNEGD R3, R5 3 and return (-Q, -ROOT) 
06 11 0040 158 BRB RETRN 
O04F 159 
OO4F 160 RETRN1: 
50 55 7D OO04F 161 MOva R5, RO 3 continue to guep ROOT and Q 
3 3 Bpae \o¢ MOVQ R3, R5 3; and return (Q, ROOT) 
005 163 RETRN: 
52 04 AC 00 0055 164 MOVL result(AP), R2 ; result address 
82 50 7b 0059 165 MOVQ RO, (R2)+ 3; real part 
62 55 7D OOSC 166 MOVva RS, (R2) : imaginary part 
04 OOSF 167 RET 
0060 168 
0060 169 
0060 170 .END 


— 
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Symbol table 66-SEP-1984 11:20:58 CMTHRTL.SRCIMTHCDSQRT.MAR; 1 (4)| 
ARG = 00000008 | 
MTHSSJACKET_HND teeeeees x 
MTHSCDABS rererere 0 
MTHSCDSQRT 00000000 RG. 01 | 
MTHSDSQRT_RS eeeneeee 0 
RESULT = 00000004 | 
RETRN 00000055 R 4 
RETRNI OOO0004F R 1 | 
¢esereeoe eee merece Sa 
} ;_Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
- (ABS. 00000000 <¢ 0.) 00 ¢ OO.) NOPIC USR CON ABS LCL — — NORD NOWRT NOVEC BYTE 
_MTHSCODE 00000060 ¢ 96.) O1¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
$onncssnenenmmaneenneand 
! ; Performance indicators ! 
Page faults CPU Time Elapsed Time 
Initialization $3 00:00:00.09 00:00:00.39 
Command processing 133 ses + SR BO De ge 
Pass 1 85 0:00:00.67 00:00:04.35 
Synbol table sort 0 00:00:00.00 00:00:00.00 
Pass 2 45 00:00:00.57 Bo 0 06 «Be 
Symbol table output ¢ 0U:00:00.02 00:00:00.03 
Psect synopsis output 00:00:00.02 BP +8 288 be 
Cross-reference output 0 hs $4 44 0:00:00.00 
Assembler run totals 299 00:00:02.0 00:00:11.92 


bytes (6 pages) of virtual memory were used to buffer the intermediate code. 
There were 10 pages of symbol table space allocated to hold 8 non-local and 2 local symbols. 
source Lines were read in Pass 1, producing 11 object records in Pass 2 
1 page of virtual memory was used to “define 1 macro. 


wee moor eno mew meno ene cnecoea $ 


Macro Library name Macros defined 
“S255SDUA2ZB:(SYSLIBISTARLET.MLB;2 0 SS” 0 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


The working set Limit was 900 pages. 
288 
MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHCDSQRT/OBJ=0BJ$:MTHCDSQRT segs <tc aadapaictie 
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